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Summary 

Infrared Earth horizon sensors in combination with a sun sensor 
have proven useful for autonomous station keeping of geosynchronous 
satellites but the complexity of a fully self-contained autonomous naviga- 
tion system for low altitude satellites has discouraged implementation of 
such a scheme. A relatively simple system which would use horizon 
crossing indicators, a sun sensor, a quartz oscillator, and a micropro- 
grammed computer is being studied. 

The sensor combination is required only to effectively measure the 
angle between the centers of the Earth and the Sun. Simulations for a 
particular orbit indicate that 2 km r. m. s. orbit determination uncertainties 
may be expected from a system with 0. °06 measurement uncertainty. A 
key finding is that knowledge of the satellite orbit plane orientation can be 
maintained to this level because of the annual motion of the Sun and the 
predictable effects of Earth oblateness. The basic system described above 
can be updated periodically by transits of the Moon through the IR horizon 
crossing indicator fields of view. The extent to which these conclusions 
may be applied to a larger class of satellite orbits is under study. 
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Introduction 

Previous autonomous navigation schemes (references 1 and 2) have 
had two characteristics which have caused them to be noncompetitive with 
normal ground navigation techniques; they tend to be low accuracy systems, 
yet they require inordinate onboard processing capability. Higher accuracy 
autonomous systems such as the Space Sextant are even more complex and 
require a large sacrifice of payload capability to perform the optical measure- 
ments and complex data reductions. To date, the most successful applications 
of the autonomous navigation concept have been for limited functions, most 
notably automatic longitude station keeping of geosynchronous satellites 
LES 6, LES 8, and LES 9 (references 3 and 4). 

The present application of interest is for a self-contained low 
accuracy ( 12 k. m.,3<7) system with minimal payload allocation requirements. 
The success of this approach hinges less on accuracy than on degree of 
autonomy and simplicity. The trap we wish to avoid is the common one of 
proposing a massive and complex system that is able to overcome all possible 
problems other than those of cost, practicality, and self-sufficiency. 

As envisioned, the completely self-contained on board navigation sys- 
tem will use one or more IR Earth horizon crossing indicators, a Sun sensor, a 
quartz oscillator, and a microprogrammed computer to deliver the 
desired overall orbit position accuracy of 12 k. m. ,3a t or better throughout 
a six month lifetime mission. Such a system has the potential to provide 
this level of self-contained autonomous navigation accuracy over very long 
mission lifetimes measured in years instead of months. It is important to 
keep in mind that the proposed system is truly autonomous in the sense 
that it is independent of other systems such as ground or orbiting radip 
beacons which are susceptible to jamming or destruction. 

To date, a 470 km circular orbit with 34° inclination has been 
studied using a special version of the FLEXSAT program. 

A brief description of this program is given in Appendix B. 
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Sensor Measurement System 

Figure 1 illustrates a conventional attitude sensor configuration 
that is well suited to perform autonomous navigation functions. The 
spinning satellite uses one or more narrow angle IR horizon crossing 
indicators and a wide angle Sun sensor. For autonomous navigation with 
horizon crossing indicators it is desirable to orient the spacecraft spin 
axis normal to the orbit plane as shown. This may be controlled by monitoring 
attitude throughout the orbital period and minimizing variations in the horizon 
scanner pulse widths by means of attitude maneuvers when required. The 
attitude measurements allow determination of the direction to the center 
of the Earth with respect to the Sun at each scan. As indicated, the hori- 
zon sensors can also detect the Moon. This opportunity will occur at least 
twice in a sidereal month. The moon observations provide an inertial refer- 
ence update that normally would require the extra complexity of a separate 
star sensor system. For the system shown in figure 1, a wide angle sun 
sensor is used to measure the times of Sun crossings through the instru- 
ment field of view and the elevation of the sun with respect to the optical 
axis of the sun sensor. The horizon and Sun transit times, along with the 
Sun elevation, yield the angle between the centers of Earth and Sun as seen 
from the satellite. 

An ambiguity exists in this measurement system, in that a rotation 
of the satellite orbit plane about the Earth-Sun line would be undetectable 
in the observations if the gravitational potential field of the Earth were that 
of a sphere rather than that of an oblate spheroid, and if the direction of 
the Sun in inertial space were fixed. The proposed system takes advantage 
of the known nature of Earth oblateness effects (see Appendix A) and of 
the orbital motion of Earth in the plane of the ecliptic. The dynamical effects of 
oblateness include regression of the nodes along the equator; the orbital motion of 

Earth defines the ecliptic plane. Periodic Moon observations remove any 
remaining ambiguity. Initial orbit knowledge at time of orbit injection 
should be sufficiently accurate (528 meters, 0. 61 m/sec ) to provide 
confidence that the ambiguity will not be a problem in practice. 
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Error Model 


For the initial studies, the sensor measurements have been simu- 
lated in the form of angular distance between the centers of Earth and Sun 
at one minute intervals during the portion of the orbit in which the Sun is 
visible to the satellite. An uncertainty of .1700 m (reference 5) was 
assumed for the uncertainty in the height of the 14-16 micron absorption 
layer of Earth's atmosphere. This translates to a horizon sensor angular 
measurement error of 

= of 042 

for a 470 km altitude orbit. 

The Sun sensor can measure angular position of the Sun to 

= of 03 
*s 

and the angular uncertainty between the optical axes of the horizon and 
Sun sensors is o 



We consider time-tag uncertainties resulting from instrumental delay and 
clock error to be similar in magnitude to <x a . The uncertainty, <j a in 
the angle between the center of Earth and center of Sim is approximately 
the r. s. s. of these errors or 

a e = { \ + \ + 2 \ ] 1/2 * °° 06 ’ 

which is the angular uncertainty used in the simulations. 

FLEXSAT was used to generate state vector covariance matrices 
based on the angular measurement uncertainties. The ballistic drag 
value, Cj^A/W, was also estimated. In addition, the Kalman filter perfor- 
mance was tested by perturbing the initial values of the estimated parameters 
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bv the amount of the a priori uncertainties. These uncertainties are 

listed in Table 1. The reference trajectory value of C d A/W was 0.037 m /kg. 

Additional filter errors were introduced by modelling an eighth 
degree, eighth order geopotential field in the numerical integration of 
the reference trajectory used to generate the simulated observations, 
whereas a second degree, zero order fit model was used. Corresponding 
covariance uncertainties were roughly approximated by adding process 
noise to the covariances in the form of acceleration uncertainties, 

<r< ^c, y, &* = ®o.2^g* 

to represent high frequency geopotential accelerations and unmodeled 
aerodynamic drag variations. The low frequency term, in contrast, 
produces accelerations of up to about lOOOfJ-g. The simple analytic dis- 
turbing function of Appendix A serves to model the accelerations very 
precisely. The velocity vector of a satellite in an inclined orbit is there- 
fore surprisingly determinable in equatorial coordinates without a stellar 
reference. 

Effect of Orbit/Sun Geometry 

Figure 2 illustrates the various possible extremes of geometry for 
a 34° inclination orbit, depending upon the time of day of launch and the 
time of year. Consider the Sun 1, Sun 2 and N axes to be in the plane 
of the drawing. The Sun 3 axis, equatorial plane, and satellite orbit plane 
are normal to the plane of the drawing. As shown, the Sun can be within 
+ 23?5 of the equator, depending upon the time of year. Sun 1 and Sun 2 
positions are extremes of solar declination. In the drawing they are placed 
normal to the satellite line of nodes so that at Sun 2 the maximum angle of the 
orbit plane to the sun line of 57°. 5 is attained. That this is unfavorable geometry 
is evident. At the limit, 90° is singular, for if the satellite attempted 
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to navigate by observing a celestial object at N, normal to the orbit plane , 
it is seen that in a circular orbit the Earth-object angle would not change as 
a function of time, to first order. 

The Sun 3 geometry is also unfavorable, as the sunline is coplanar 
with the satellite orbit. Since this particular configuration again places 
the satellite line of nodes on the ecliptic line of nodes, the orbital inclina- 
tion of the satellite is not directly observable. 

Position Uncertainties 

Table 2 contains the peak remaining radial, intrack, and crosstrack 
position standard deviations for each of the three extreme Sun orientations 
after nine simulated orbital revolutions of the satellite, using FLEXSAT 
covariances. An advantage of a recursive real time filter is that the 
customary prediction errors are limited to data gaps, which in this case 
are somewhat less than half of each orbit revolution. 

The largest crosstrack errors are associated with Sun 1 orientation, 
with Sun 3 a close contender. Figures 3 and 4 plot the time history of 
these covariance -derived uncertainties as a function of time from injection. 

It is seen that the orbit solutions are stable but not overly convergent. 
Simulations with a spherical Earth model produce crosstrack uncertainties 
that increase with time, as expected, in the presence of the 80jUg acceleration 
noise that simulates unmodeled high frequency geopotential and drag terms. 
This results from the ambiguity in the orientation of the orbit plane that 
would exist except for the measurable presence of the disturbing function. 

The largest position uncertainty was found to be the intrack standard 
deviation in the Sun 2 configuration, which also produces the largest radial 
uncertainty (see Table 2). The time history from injection of these errors 
is plotted in figures 5 and 6. As expected, these exhibit more convergent 
behavior than do the crosstrack uncertainties, which more closely reflect 
orbit plane orientation errors. However, it should be noted that the highest 
intrack errors are initially large and do not converge to the extent of 
recovering a priori knowledge. The Earth horizon measurement errors, of 
course, map directly into intrack orbit errors. 
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Effect of Injection Knowledge 

The previously described cases used the somewhat conservative 
orbit injection knowledge uncertainties and errors in Table 1 of 


7f - H 

^x,y,z a x,y,z 

• • S ft* • • 

<7x,y,z "x,y, z 


3048 m 
3. 048 m/sec. 


The crosstrack errors for a typical nominal case (Sun at first point of Aries, 

Q = 90°) are plotted from injection through nine orbit revolutions in 
Sat 

Figure 7, 

To verify the dependence of orbit plane orientation knowledge on 
injection knowledge and to test the capability of retaining this knowledge, a 
similar case with more realistic injection knowledge and errors was run 


using 


<7 x,y, z, 4 x,y,z 

a, . . = 2. . . 

k, y , & , x, y, z 


528 m 

0. 61 m /sec. 


Radial uncertainty was reduced from about 600 m (nominal case) to 400 m 
on the ninth orbit revolution, while intrack uncertainty was reduced from 
2700m (nominal case)to 1900m on the ninth revolution. The important 
crosstrack uncertainty is plotted in Figure 8. It is seen that the injection 
knowledge of 528 m is retained through the ninth orbit revolution and even 
improved slightly between the first and ninth revolutions. This is certainly 
encouraging in light of the importance of minimizing orbit plane orientation 
errors. As expected, however, the solution displays slightly divergent 
characteristics. In time the errors might be expected to grow to the size 
of those in Figure 7. 

Filter Errors 

The largest filter estimate difference from the "truth" state vector 
or from the "truth" ballistic drag value is less than 3<T, where the value 
of <7 is obtained from the covariance matrix associated with the particular 
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estimate. Assuming a Gaussian error distribution, one would normally 
expect to see an occasional 3ff estimate. The great majority of estimates 
are less than Iff from the "truth" model. The actual estimates are 
accordingly better than Table 2 and Figures 3-8 indicate. In these cases 
the conservative process noise of 80. 2/yg served to maintain the filter 
covariance matrix at a reasonably high level. The fact that there was 
convergence and that the actual errors showed reasonable conformity 
with the sigmas from the covariance analysis indicates that the filter 
covariances are realistic. Since in theory they represent an infinite sample 
of Monte Carlo trials they are the numbers tabulated and plotted in this 
paper. 

Moon Observations 

Figure 1 illustrates how the horizon crossing indicator will, in 
general, view two portions of the Moon's orbit (the second view area is 
on the opposite sides of the satellite and Moon orbits). When the Moon 
enters these view areas, once every sidereal month for each portion, the 
1^-16 micron bandwidth horizon sensor will detect the Moon for several 
satellite orbit revolutions on each occasion. The exact length of viewing 
time depends upon the horizon sensor field of view and the inclination of 
the satellite orbit plane to the orbit plane of the Moon. 

The Moon observations can be used to periodically update the orbit 
knowledge with independent observations. These observations fix the 
satellite state in inertial space in a direct manner. If two horizon crossing 
indicators are used in order to scan both north and south of the orbit plane, 
then two additional Moon viewing periods are available in each sidereal 
month. This system would seem to be superior to a system using 
only one horizon crossing indicator in any event, when the attitude 
determination problem is examined. 

The principal value of Moon observations is to provide periodic 
recovery capability in the event that orbit knowledge is lost or degraded 
owing to larger than expected injection errors, degraded sun sensor per- 
formance, transient data stream/clock/microprocessor failures or un- 
expectedly large perturbations to the satellite orbit. A very compact, 
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pre -calculated lunar ephemeris would suffice because of the infrequency of 
Moon observations. With two horizon sensors, the Moon would nominally 
be observable for four or five orbit revolutions per week. To account for 
the large IR radiation differences between the illuminated hemisphere of 
the Moon and the dark side, a lunar phase -dependent model of the asym- 
metric sensor response would be a necessary part of the pre -calculated 
ephemeris. Figure 9 shows the essential elements of the navigation system, 
including the Moon data capability. The dashed lines indicate that the Sun 
elevation measurements are optional for attitude control, but may prove 
useful. 


Clock Errors 

All horizon, Sun, and Moon observations must be time-tagged by the 

9 

onboard oscillator. A typical quartz oscillator is stable to one part in 10 , 
or 30 msec /year. To make use of 500 m injection accuracy we desire 
clock errors no larger than 


. 500 jc 5640 

m ~ 6848252 x 2 7T 


» 66 msec 


during the intervals between Moon observations, which is clearly not a 
problem. In the above example the orbital period is 5640 seconds and the 
orbit semimajor axis is 6848252 m. 


Onboard Computer Requirements 

The products of the autonomous system diagrammed in Figure 9 are 
the satellite ephemeris at bottom center and the attitude control function at 
upper right of the chart. The recursive orbit filter and attitude computations, 
ephemeris evaluation, and information management throughout the system could 
be performed by a microprogrammed I/O and central processor system. 

The requirements are currently being studied, but it is estimated that a 32k word 
memory and 16 bit fixed word length should be adequate. 
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Conclusions 

An autonomous navigation system such as the one diagrammed in 
Figure 9 would appear to be capable of delivering accuracies normally 
associated with horizon sensors in conjunction with a stellar attitude 
system (reference 7). It is felt that the extra complication of a star 
sensor may be unwarranted considering the relatively good performance 
of a horizon sensor and sun sensor system. To fully assess the value of 
such a system., however, it is important to study the particular orbital 
characteristics of the intended mission. For example, simulations indicate 
that some high inclination missions may be a poor choice or would at 
least require further study. Accuracy will also be dependent to some 
extent on orbit altitude. 
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Appendix A 


The disturbing function of an equatorial bulge is 

D = J-, (3 sin 2 <J -1) 

2r^ 2 

where U is GM ?apfch , 

r is the instantaneous radius vector, 

J 2 is the second degree Legendre polynomial coefficient 
for Earth, 

6 is the instantaneous declination of the satellite. 


The secular perturbations are then 

1 r zn 

- TWJ 0 D d M 

where M is the mean anomaly. 

In terms of orbital elements a, e, i, M 0 ,the principal secular 

effect is a regression of the nodes along the equator, 

3n J,, 

cos i dt. 


dQ s = - 


o 2/- 2,2 

2a (l-e ) 


where n is the mean motion. 

Depending upon whether orbital inclination is less than or greater than 
i = arc sin (2 Jy/~5) - 63. %3 , the line of apsides will secularly advance 

or regress according to 


3n J„ 


d Id = 
s 


2 2 ( \ sin^i -2) dt. 


2a (l-e ) 

Secular changes in the orbital period are also a function of a, e, and i as 
the mean anomaly changes by 


3 J, 


d M = n d t [ 1 - 
s L 


2a ( l-e ) 


,3.2. 

2.3/2 ^ 2 Sln 1 " 



APPENDIX B 


SIMULATION SOFTWARE 


I program-. FLEXSAT 
I author: H. Hendrickson 

I compiles on: CDC 7600, CDC 176 FORTRAN compiler 

• trajectory: Runge-Kutta fourth order integrator 

• filter: Kalman 

• partial derivatives*, finite difference 

- f (p^, . . . . , p j + 4Pj, . . . P|^)({.) “ f (Pji • • • • * P|(\tj) 

at the ith observation, 0j and jth estimated parameter, 
Pj at time tj for the data equation 

( Pl ,Pk )+ € 

having k parameters, e observation noise. 

• Additive parameter noise model: 

White acceleration noise, 

.7 = O I (t. -t. ) 1/2 . 

q. i+l i 

r r 


A 


mil 


III 


ONBOARD CALCULATIONS FOR HORIZON/SUN SENSOR AUTONOMOUS 
ATTITUDE AND NAVIGATION SYSTEM 


I Polynomial evaluation, one second intervals for observation angles, 

3 Annual Terms e (a, e, M ) 

t o Earth 

8 Monthly Terms 0 (a, e, i, Q,u, M , 

(use different a Q t * 0 

a u each month) aq,au) 

Moon 


3 Secular terms 0 (aq, aoj, aN\) 

(Satellite orbit) * J 2 0^2 2 

I Sequential filter recursive estimation at each observation time 

2 attitude angles ( a, 6 spin axis) 

3 satellite orbit parameters 

(a, e, i, q, cj, C A/W) and ballistic drag 

D Satellite 

(assume M 0 known at injection) 

• Satellite ephemeris calculation 

geocentric position 

as a function of time (x, y, z) 

at one minute intervals * 



Simulation a priori Uncertainties/! nitial Perturbations 
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